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A STUDY OF BLAST LOADING 


TRANSMITTED TO BUILDING FRAMES 


1. INTRODUCTION 
A. Introductory Statement 

The loading which a structure experiences due to the pressure waves accom- 
panying blasts from high explosive and atomic bombs has been the subject of 
much study and research since the use of such bombs vas begun in the last var. 
The advent of more poverful atomic bombs and of the "super" bombs has increased 
even more the urgency of discovering vith greater exactness the nature of blast 
loads on structures. 

A recent study in the field of blast loading ol” has concluded that a gener- 
alized triangular pulse is a good approximation to the actual transient loads 
which the front wall of a structure experiences. However depending on the 
action of the wall covering, the shape of the load-time curves that reach the 
building frame itself may be quite different from the pulse which impinges on 
the wall panels. It is the purpose of this thesis to study the load-time rela- 
tionship for the load transmitted by the wall covering to the supporting frame. 

The problem is necessarily complex and is approached by making use of cer- 
tain simplifying assumptions. The most important are: 

(1) Given a typical bay of a one-story frame-type buildíng, consider that 
the wall covering of the bay consists of a "nest" of closely-spaced 
beams of unit width which are simply supported and act independently 
of one another. (Figure 1) These beams are used in this study to 
approximate the action of a slab covering. 


DI M. Nevmark, "An Engineering Approach to Blast Resistant Design’, 
Proceedings ASCE, Vol. 79, Separate No. 306, October 1953. 


ze 





(2) Assume that the pressure loading which impinges on the wall panels 
can be represented by a triangular pulse rot (Figure 2) 

(3) Assume ideal elasto-plastic action of the beams. (Figure 3) 

(4) Assume that the frame members supporting the wall covering remain 
fixed in space, at least during the interval for which this study 
is made, 

It is doubtful that movement of the frame would have much effect on the 
loads reaching the frame members. A deflection of the frame away from the 
blast source would result in a slight decrease in the beam end reactions; 
but, since the period of vibration of the supports is much larger than that 
of the wall panels, the effect should be negligible except possibly for load 
pulses of very long duration. 

If the beams which approximate the wall covering were perfectly rigid and 
massless, the shape of the load transmitted to the frame would be identical 
to that of the blast pressure pulse (Figure 2). However in our case, the beams 
have mass and also have elasto-plastic characteristics; therefore, the load 
which reaches the supporting members is much altered. 

In order to determine the frame loading for the given assumptions, an 
elastic Keyat and a simple plastic analysis are used. The elastíc analy- 
sis is derived from the familiar dynamic beam equation (1) and the plastic 
one comes from a simplified consideration of beam action when a plastic hinge 
forms at the center. From these two analyses expressions for the dynamíc end 
reaction of the beams on the supports are derived. 

Frame load charts giving the variation of dynamic end reaction (frame load 
per unit height) with respect to tíme are drawn for different beam resistance 


om M. Newmark, "An Engineering Approach to Blast Resistant Design", 
Proceedings ASCE, Vol. 79, Separate No. 306, October 1953. 
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characteristics and different load pulse durations. The frame loads are com- 
pared graphically with the load which the beam receives as a result of the 
blast pressure wave. Load-time curves for the supporting frame are shown 
only for the duration of the loading pulse; however, a method is indicated 
by which the relation may be found for times beyond the end of the pulse. 

As brought out in a subsequent section of this study, the elastic reactíon 
expression involves the sum of a trigonometric series. For certain durations 
of loading pulses, the curves for one term, three terms and six terms are 
drawn for comparison. By means of this graphical comparison of a few repre- 
sentative pulse durations, it is demonstrated that the use of only the first 
term gives results which are fairly accurate for long pulses but not very ac- 
curate for short pulses. For all the remaining pulse durations, only the 
first mode (one term) is used although this is only an approximation to the 
exact relation which theoretically would require an infinite number of terms. 
No attempt is made by analytical means to evaluate the accuracy obtainable for 
various pulse durations by using only the first mode. Such an evaluation 
could be made the subject of further study. The advantages of using the first 
mode approximation in the construction of the frame load charts is apparent, 
however, and such use is intended. 

B. Notation 

Terms used in the following sections are defined as they appear. They 
are assembled here for convenience in reference. Dimensions of the quantities 
are indicated by the symbols F, M, L, T in parentheses meaning respectively 
force, mass, length and time. 

a? 2 EL/n (L /T”) 

~ = angular displacement 


9 = angular velocity (1/T) 





angular acceleration E 1/7”) 

impulse of load = püt (FT/L) 

time distribution factor 

SE dr.) 

length of beam (wall panel strip) (L) 

mass per unit length of beam (M) 

dynemic bending moment at center of beam (FL) 
fully plastic moment of beem (FL) 

maximum static center moment = P 2 2/8 (FL) 
ultimate moment - carrying capacity of beam (FL) 
instantaneous value of applied pressure load (F/L) 
7 Së = space distribution factor 

maximum value of applied pressure load (F/L) 
dynamic end reaction (F) 

maximum static end reaction = Pt /2 (F) 

time (T) 

duration of load pulse (T) 

time of beginning of second elastic phase (T) 
time of beginning of plastic phase (T) 

5 dü 

t-t, (T) 

natural periods of vibration of beam (T) 
fundamental natural period of vibration of beam (T) 
a function of time alone 

shear (F) 


deflection of beam (L) 





v m linear velocity of beam (L/T) 


w = distance from center of mass to axis of rotation (L) 

v = acceleration of mass center (1/Тб) 

Wo = effective center deflection of beam at onset of second 
€ 


elastic phase (L) 





w, = maximum center deflection of beam which occurs at t = t. (L) 
P: = center deflection of beam at t = ER (L) 

X. = characteristic function for beam T 

w, = natural (circular) frequencies of beam = T VE (1/T) 
W, = fundamental natural frequency of beam (1/T) 


II. INITIAL ELASTIC BEHAVIOR OF WALL PANELS 
A. Statement of the Problem 
1. The Dynamic Beam Equation cə) 
When a beam resists moving loads, the governing differential equation of 
motion for the beam is, according to the ordinary Bernoulli-Euler theory of 


flexure, 
= 


(1) 





ox öt m = 


where ER = EI/m; w = deflection measured normal to the beam axis; m = mass per 
unit length; q = external applied unit load; x = distance measured from origin 
at left end of beam; and t = time. The external applied load is defined by a 


a 
space distribution factor, p(x) « 9 P, X. and a time distribution factor, p(t). 


n 


see for example, S. Timoshenko, "Vibration Problems in Engineering”, D. Van 
Nostrand Co., Inc., New York, N. Y., Second Edition, July 1937, p- 332. 





Figure h is a simplified sketch of the loading on a typical beam. For the 


- 


given triangular loading pulse, p(x) = l since the beam is assumed to be uni- 
formly loaded, and Ø(t) = Р (1 = t/t,). The characteristic function, X, for 


a simply supported beam is sin (n x/.£ ). The constants, Pi are defined as: 


Can x. ds 
P. = 7 (2) 
2 
L vs 


Substituting in (2) and performing the integrations leads to: 


o A 


I! 


h 
— for n 
n Ол 


"Ú 
N 


A 00 A en 


The solution to (1) is taken in the form 


where X, = sin (nrx/£ ), a function of x alone, and where pa is a function of 
time alone. The substitution of (3) in (1) vith the subsequent separation of 


variables leads to the folloving differential equations: 


a 
Да (s) 
a 
L 
1 0 _ Py Plt) (5) 
n n n m 


The solutions of (4) and (5) are respectively: 





Y 
X = sin T (always true for simply- (6) 


supported beam) 


T zx A cos Wt +B sin wt 
n n n n n 


P t k 
E. sin ant f P(t) cos w tdt - cos dit Ö(t) sin w tat (7) 
n O o 


Eguation (7) is the solution to the single-degree of freedom Equation (63.90 
The first tvo terms in (7) are the free vibration terms and the last tvo the 
forced vibration terms. For initial conditions of zero displacement and zero 


velocity, the free vibration terms disappear since 


where v (x) and vo (x) are the displacement and velocity respectively at zero 


time. 


Js. Timoshenko, "Vibration Problems in Engineering”, D. Van Nostrand Co., Inc., 
New York, N. Y,, Second Edition, July 1937, p. 348. 





8 
Substituting Equations (6) and (7) into the expression for deflection (3), 


yields the following: 


P 


n 
Ku —— Sin 
mw 

n 


t t 
Dax 
7 sin wt D(t) cos w tdt - cos w t ö(t) sin u tdt (8) 


o O 


Substituting the values of P_ and Ö(t) in (8) and integrating, results in the 


following expression for w: 





1 m nx t sin әј 
w = aa o ` ae (1 - cos wt)- += 1- q 
non 1 n 
n 
2 2 2 
Noting that w for a simply-supported beam = (n n/ LT) a, where a =,/ El/m: 
up 23 sin o t 
m 1 nT x t n 
wa — “5 sin = (1 - cos wt) ə € - (9) 
"EI n 1: n 
"ијә => 


B. Reaction and Moment Expressions 


1. Dynamic End Reaction 
The shear at any point x along the beam is given by: 


3 
v - - El Sİ 
әх? 


Performing the differentiation on (9) and evaluating the expression at the 


end of the beam (x = O), the dynamic end reaction is: 








| iz 1 t sin wt 
Cher SS A EL ee 
TT n 1 Ən 


n=l,3,5--- 


If — the maximum value of the external load, were applied as a static load, 


the maximum statíc end reaction would be: 


For ease in plotting, the reaction expression is put into dimensionless form 


2 2 
by noting that ә an/T = 2n 2 since Ta = GIG s 





2 
sin 2n“ t/T 
To = > 1 - cos nn” = - 1/(t,/T,) ° t/T, 1 - 5 1 10) 
Rg n n 1 aun t/T, 
n=1,3,5--- 


2. Magnitude of Reaction for Long Pulses 


As the duration of the load pulse goes to infinity, the coefficient of the 


second term in (10) approaches zero and the reaction approaches: 


Бб 1 . 

lim == = = l - cos mn = 

Rs ef ne 1, 
n=1,3,5--- 


Be => ei 


The maximm reaction then approaches: 
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oo * “ 
8 2 16 it 1 -X 
lim 5 — = lim > -5 , where lim “ee bez 6 (x) 
D n TI n n 
n-odd n=odd n=odd 


The limit of the summation exists as shown. It is given by the Riemann Zeta 
Function, the values of which are Ein. °°? 


ә. Ele) = "5/6 = 16645; or: 


кә 
2 
tin Ë e 5 5-8 (1 - 1/4) > = 2.00 
TI n Y 
> e n=odd 


Then for pulses of duration approaching infinity, the dynamic reaction ap- 
proaches a value twice the maximum static reaction. This result is to be 
expected from energy considerations. 


3. Magnitude of Reaction for Short Pulses and Constant Impulse 





For short pulses of duration approaching zero, the reaction should also 
approach a limit. This limit can be found from a consideration of the reac- 
tion expression which holds for time beyond the end of the pulse as the pulse 
duration becomes very small and total impulse remains constant. Figure 5 


shows the new time origin at t = Eu: 


As before, the deflection, w = > Х, тр" However for t > ty the forced 


n 
vibration terms in t, drop out, leaving: 
x x 
=A_ Cos wt +B sinut 
na n n n 


T 
n 


> Y 
m Jahnke and F. Emde, "Tables of Functions with Formulae and Curves’, 
Dover Publications, New York, N. Y., Fourth Edition, 1943. 





7 
The initial conditions of displacement and velocity for t - 0 can be found 


from Equation (9). Setting t - t, in (9): 





LP 48 sino t 
w (x) 2 S SC Sin T E - COS ü у 
t “El n “o 1 S 
n=odd 








1: 
Lp L" 1 dalı cos Mot? 1 
v (x) = — sin | sinw t., + = 
: WEI -. L — E = 
n= 
Values of An and BA are then: 
7 v (x) X dx 
a 5 e nux 
A = ——— where X = sin — 
n 4, n £ 
2 
f <a 
o n 
up Di sin ut. - uk co 
2:50. az nr an n l nu 
a "?ЕТ = SN 
n=odd 











dı 
h 
y up Ê 8 wt, sin wt, + cos wt, - 1 
= rl A EE 
D PE n? e 
n=odd 
x 
Then for t = t - 5): 
- 
4 
up £ Sinwt. - wt. cos wt 
w(x) = m - Sin 7 cos ut” — 1221 
n El n n 1 
n= 
E yt sin ә 6) + cos Wu - 1 
n 1 


As before, find the shear and evaluate at x = O in order to get the dynamic end 


reaction: 





ir. 2 i „[ein ot, - wt, cos ut, 
Rp 2 2 m W 
T n nl 
n=odd 
Wt. gin ú t. +O cos u t. - 1 
+ sin ort ñ l no nt (11) 
n w t 
n 1 


For any pulse duration, the impulse of the load is: 


t1 
I = d pdt = 1/2 Paty 





Now let tl —> O keeping the total impulse constant. 


l 
rapidly, and in terms of I is given by: 





Replacing P| in (11) by its value in terms of I and taking t 

mation: 

sinWt, - wt 
n n 


2 
N 


“o 


cos uwt l 
n 


1 





u) 
к 


sin Wt. + cos Wt, - 
n 1 n 

2 

u) 

d 


+ sin Q t 1 
n 


Now let t 


1 approach zero: 








Noting again that u = an/T and replacing I by (1/2) Pit 


oo 
x 
8t sin u t 
_ 1 5:5. 
"n | nT 
n 


n=odd 


i 


R 
Lim == 
5 


t О 


1 


il 


Take the first mode by letting n 


и 


13 


Ер then increases 


1 inside the sum- 


cos a 





1) 


lim 
t, — o 


= | O 
ct 

c+ 
` 


AS 


1 
E sin 21 m (12) 


The foregoing derivation has shovn that for the conditions of very short 
pulse duration (say vhen b < 1/3) and constant impulse, the maximum elastic 
reaction approaches a value equal to (8/n) t,/T, - A plot of reaction for 
t > t, in which the load pulse satisfies the two conditions should verify (12). 


h. Dynamic Center Moment and Beam Resistance Parameter 
The moment at any point x along the beam is given by: 


2 
w 


Q/ 


Mos - El 


| 


e 


Performing the differentiation on (9), evaluating at x = //2 and noting that 
2 
Ef /8 is the maximum static center moment if P were applied as a static 


load: 


n-1 
M, a) è 26 
M” CR 1 - cos enn m 
n=odd 
sin anne ə 
farm) um (1 — (13) 


A plot of мм, versus the time parameter t/T, furnishes a convenient means 
of defining the limits of elastic and plastic beam action for different pulse 
durations. As in the reaction plots, the first mode (n = 1) in (13) is used 


as a good approximation. Given a curve of M/M versus t/T, , the fully plastic 





15 
moment of the beam cross-section may be considered to have been reached for 
different values of Мр. When a given MYM. is said to equal M.M, the time 
corresponding to the particular value of Мм. is the time of hinge formation, 
— The beam resistance parameter, M /Mç may be thought of as the ratio of 


the yield value of the beam's moment resistance to the maximum moment caused 


by the static application of ba: 


III PLASTIC BEHAVIOR OF WALL PANELS 
A. Statement of the Problem 
1. Assumptions in the Simplified Plastic Analysis 

In order to derive a fairly simple expression for beam reaction when the 

beam ceases tc behave elastically, the following assumptions are made: 

a. A plastic hinge forms at the center of the beam when the dynamic center 
moment , Му reaches a value equal to the fully plastic moment of the sec- 
tion, Мр. This time is designated E 

b. Neglect the elastic bending of the bean. 

c. Assume that deflections are small. 

d. Assume that Му remains equal to Mp until the beam velocity becomes 
equal to zero. At this point, Му becomes less than Mp» the center 
hinge disappears, and the beam behaves elastically again. 

Figure Ó represents the idealized plastic behavior of the beem conforming 

to the above assumptions. 

Figure 7 is a free body diagram of half of the beam from which equations 


of motion may be written as follows: 


= SC E S.S 
| F. = Mw where w = f 9; Wwe © 
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fonol Ea Qu) 
> 2 3 
A big where I E nl (4) .. 
Sg 
cəkən. üs) 


B. Reaction Expression 


From Equation (15) solve for the angular acceleration: 
Substitute (16) into Equation (14) and solve for R: 
R = SW + ” vhere p function of time. (17) 
In Equation (17), p is the external load acting on the beam at time t" measured 


after the formation of the center hinge. Figure 8 represents the time origins 


for the different phases of beam action vhich occur during the load pulse. 


t' + t 
p= p(t') =P | 1 - | —— 
m t 


and t' = t - ty and noting again that P, (£/2 is the maximum static reaction: 


With 
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Ер 1 J EN AE 24M, 
= 1 - 1/(t./T — €, (18) 
Ro D 771 бу, pL? 


Differentiation of (18) with respect to the time parameter, Gas gives the 
result that the slope, S, of the plastic reaction curve depends only on the 


pulse duration, or: 


S=- r tan. = -В if Ke 1 (t,7,) (19) 


C. Duration of Plastic Behavior 


1. Expression for to 


As stated in the assumptions for the plastic analysis, the center hinge 
disappears when the beam velocity becomes equal to zero. The time, to? at 
which the beam stops moving down is therefore the end of the plastic phase 
and the beginning of the second elastic phase. 

In order to find to an expression for the angular velocity of the beam 
can be derived which holds during the plastic phase BE SE ta). Equating 
this angular velocity to zero and solving for the time will yield the equa- 
tion for tə 

Differentiating Equation (9) with respect to time and evaluating at 
x = 4/2 gives the linear velocity of the beam center, y, during the initial 
elastic phase. For small deflections, the average angular velocity of half 
of the beam would be = 2/2 “. Performing these operations and Setting 
“ölkə Sr will give the equation for the initiel angular velocity at t' = O 


(t = 5). (Figure 8) Thus: 
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° E 
TI 


Bo 52 
8P 
gr = u w sin ut - = (1 - cog o t ) (20) 
El n P 1 na 
odd 


The required angular velocity may nov be found by integrating (16) vhich 
was derived from the equations of motion for the plastic phase. 


Noting that p in (16) = 


"at 
p(t') = P, (CZ 


t' 
s - m A NÉE 
m Í l PL 


Integrating and evaluating the constant of integration by using (20): 


: бкК.. . ı ft: Mp 


LA 


XV. SS 
Em L (-1) 1 ( 
+ w sin wt_ - =— [1 - cos u t 
“dı ) n? n n p ty np 
n=odd 


4 
Simplify by letting n = 1 (first mode) and note that ә = ""El/m2 : 





1, 








P 
Ds | a E) `" 
t 2 p M 
mL 1 S 
+ = (1 - cos ut ) (21) 
3nuy lp Wt, l p 


When (21) is equated to zero, the resulting equation is a quadratic in t' 
which can be solved by the quadratic formula giving the following result when 


put in dimensionless form: 





3 l on E 


Ae 


IV. SECOND PHASE OF ELASTIC WALL PANEL BEHAVIOR 
A. Statement of the Problem 
1. Application of Dynamic Beam Equation 
Equation (22) defines the time at vhich the beam again behaves elastically. 


fhe derivation of this second elastic reaction expression follows the same 
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method as used on the initial one; however, one free vibration term is added. 
When t" = O (t = to)» the initial velocity, e. of the beam is zero, but the 
initial displacement, v has a finite value, An examination of the formulae 
for A, and B, (Section II) shows that the first free vibration term in Equation 
(7) has a value other than zero. 


For time measured from t, (t" = t - te), the elastic deflection is still 


of the form: 


w"(x) = 1 Х ту 


n 
Where X, = sin nand and 


F 
7 =A cés wut" — sin ú vr P(t") cos w t'at" 
n n n mo, n J, n 


es 
- cos | ö(t") sin w t'at" (23) 
o 


where P_ = 4/nu when n = odd. (See Equation (2).) 


The constant Ад is defined as 


1 

УД w"(x) X, dx 

— s +f w_"(x) X, dx (24) 
N 


2, 
o 





eL 
Now let wo (x) = v. sin 7 where v. is the effective center deflection at 
the beginning of az second elastic See, At first glance it appears that a 
good approximation of wo (x) should be given by 


wo (x) =w sin Le 


Co 7 


иһеге "Ze, is the center deflection of the beam in the plastic phase at the 
instant the beam reaches its maximm deflection in the positive direction (in 
direction of applied load) and begins to rebound elastically. However, it is 
apparent that the maximum amplitude of the elastic center deflection is much 
less than the corresponding maximum amplitude of the plastic center deflection. 
Figure 9 illustrates roughly this difference. The solid curve gives an idea 

of the actual variation of reaction with center deflection throughout the two 
elastic phases and the single plastic phase. The dotted curve is an approxima- 
tion to the actual behavior which suggests a simple means of obtaining a close 


estimate of v. . As shown in the figure, the maximum plastic center deflec- 


e 
tion, Wo is much greater than 2: and m , the center deflection at the time 
m e 
of hinge formation Boch is approximately equal tow, . 


e 
The final expression and best approximation to be used for w (x) is thus: 


wo (x) =w sin €” 


y £ 


The value of v, may be found from Equation (9) by setting x - //2 and 


t = Se: The result is: 


n-1 
SÉ 


sin w E 


, e 
hp / t 
m -l n 
W = Cu (1 - cos ut )--E 1 - —_ (25) 
У wiet | n "P 2 ( va p 
n=0dd 
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nx 
Now ín X = sin a (24), let n= 1. A, then reduces to: 


L 
2 2 TX 
El EE тк. 


The final expression for A is thus identical to Wy and is given by (25). 


Go back to the expression for R. Equation (23), and substitute the value 


of Alt") 


27—:42:-.— 


After carrying out the integrations in (23), combining it vith X, end simpli- 





fying, the expression for elastic deflection for Ки кисә ty is: 
00 
иј? 1 nT x 
w"(x) = — sin Fr cos w t" (1 - coğbu 6.) 
"ЕТ n? 2 SE 
n=odd 
t sin wt t. sin ut” br t. 
ee AAA H (00007 AAA ES) 
t u t t w t t 
1 n 1 n 1 1 


B. Reaction Expression 
Having derived Equation (26) for deflection, the reaction expression for 
the second elastic phase may be found as before by evaluating the shear at 
x = O. In dimensionless form for plotting, and using the first mode only, the 


reaction equation for t. € t < t), is: 





e3 


| t Ë schuet, Wm 
R, 8 = P 1 P al 
— - — ( cos 21 — 1 - cos 27 - . 1 - 
Ro ne Т) Т, Enn Т) on t Tn 
t sin 2n t"/T tu 
Birk) Pr) |—— (27) 
r” 1 51 1 1 


V. RESULTS AND DISCUSSION 
A. Frame Load Charts 

1. Construction of Charts 

In the previous sections of this paper, expressions for the beam end re- 
actions have been derived for the elastic and plastic phases of beam behavior. 
The reaction is expressed as a ratio of the dynamic end reaction to the maxi- 
mum static end reaction, (RR); the time parameter is written as the ratio 
of the elapsed time to the fundamental natural period of the beam ((/1))) and 
the beam resistance parameter is expressed as the ratio of the fully plastic 
moment of the section to the maximum static center moment (M,/M,.). For con- 
venience, the following table lists the various phases of beam action, the 


time range of each phase (Figure 8), and the governing reaction equation: 


Phase Time Range Equation 
Initial Elastic ge te t. (10) 
Plastic CELTE (18) 

p ë 
Second Elastic t,<t<t, (27) 


Figures 23, 27, and 29 contain plots of all three phases. The remaining 


reaction curves include only the initial elastic and the plastic phases. On 
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the same graph with each reaction plot is shown the external load which im- 
pinges on half the beam. Thus at t = O, the load on half the beam is p_ 4/2; 
the maximum static end reaction is also P £ /2; therefore the initial ratio 
of the half-beam load to the maximum static reaction is equal to one. The 
external load line thus starts at unity at t = O, and drops in a straight 
line to zero at t = t): This is, then, a graphical representation of the load 
which arrives on half of the beam, end what this half beam transmits on to 
the frame as end reaction. 

The moment curves are used to find the value of SS for a given pulse 


duration t., and resistance parameter (M,/M.). (Figure 16) The duration 


P 
of the initial elastic phase is thus determined graphically from the plot of 
M/M. versus t/T, by allowing M, to equal M, as M, changes over a range of 
values. This brings in the shape and physical characteristics of the beam 
since the fully plastic moment varies with the shape and material. 

The duration of the plastic phase is fixed by Equation (22). The loci of 
these times as found by solving (22) form a smocth curve and are plotted on the 
reaction graphs. It will be noted that for a given pulse duration, the loci of 
Equation (22) indicate the expected results; that is, that beams of low resist- 
ance (small M/M..) have long plastic phases and vice versa. 

In the construction of the reaction curves, it vas found that a discon- 
tinuity always existed between the elastic and plastic phases. This "Jump" 
in reaction value is due to the relatively rough assumptions of the plastic 
analysis used, the main error being the neglect of elastic bending of the beam. 
Since the source of the discontinuity was known, the reaction values were cor- 
rected by the amount of the initial reaction error. 


Plots of the second elastic phase were made only for a few representative 


values of tn and MAM. as follows: 
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Ef M/M, Fig. 








a. As an example of a short pulse 0.8 1.0 23 
b. As examples of long pulses 520 1.8 27 
20 1.0 27 
10.0 1.4 29 
10.0 1.0 29 

(6) 


2. Comperison of Fírst Mode and More Exact Reaction Curves 


As stated in the íntroductory statement, no attempt is made in this study 
to find out by analytical procedures the accuracy with which the first mode 
approximation represents the exact elastic reaction curve. However, for the 
initial elastic phase, Equation (10) was expanded and data was collected to 
plot curves for one term (n = 1), three terms (n = 1, 3 and 5), and 6 terms 
(n= 1, 3, 5, 7, 9 and 11) in order that a graphical comparison could be made. 

Figure 11 shows the results of the plot for a very short pulse. It is 
Obvious that the first mode approximation is not good for load pulses of such 
short duration. Figure 35 illustrates the very great improvement in the ac- 
curacy of the first mode approximation for pulses of very long duration. It 
would seem reasonable to expect then, that as the load pulse duration increases 
from a small to a large value, the accuracy of the first mode should show a 
gradual improvement. Figure 19 furnishes evidence of this trend since it shows 
the first mode curve to be a better approximation than it is for GC = 0.1 


but not as good as for t,/T, = 100.0. 


| mis subsection refers to the uncompensated first mode reaction curves as 
derived from Equation (10) by setting n = 1. These curves are indicated in 
Figures 11, 19 and 35. As brought out in the following subsection, all other 
reaction plots are constructed using the first mode multiplied by a compen- 
sating factor. 
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One situation in which the first mode curve gives a false impression is 
encountered at very small values to t/T, and M/M, - An example is Figure 11. 
Here the reaction in the plastic phase for M/M, = .05 appears to become neg- 
ative before the end of the pulse. It is apparent, however, that negative 
values would not be indicated if the plastic phase plot were started from the 
more exact initial elastic curve (for 3 or 6 terms). An examination of pulses 
of longer duration will indicate this same ambiguity at very small valves of 
x... 
3. Method Used to Make First Mode Approximation More Conservative 

Figures 11, 19, end 39 show that the first mode reactions plotted from 
Equation (10) generally furnish values lower than those indicated by the more 
exact curves. In order to retain the advantages of the first mode approxima- 
tion but at the same time to obtain more conservative end reaction magnitudes, 
the reaction expressions for the two elastic phases (Equations (10) and (27) 
with n = 1) were multiplied by a factor of a /8. This step also increased the 
plastic reaction values by the same factor since the plastic phase merges with 
the two elastic phases. The reasoning behind the use of this factor follows. 

The use of the first mode curves for elastic reaction neglects the contri- 
bution of the higher modes in (10) and (27) to the end reaction. Such a sin- 
plifying approximation may be considered as being equivalent to using an in- 
stantaneous beam loading represented by the first term in the Fourier Series 
expansion of a uniform load. This first term in the Fourier expansion is merely 
a Simple sine curve distribution of load. Therefore, the actual pressure load 
which at any instant is uniformly distributed, is replaced by a sine curve load- 
ing the total magnitude of which is very close to that of the actual load. 
Figure 36 illustrates the shapes of the actual and approximately -equivalent 


loadings. 
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The ratio of the end reaction calculated from the correct loading to that 


derived from the equivalent sine curve loading is found as follows: (Figure 36) 


Hü = 
1 h 1X hp(t) 2 
R(2) = 5 | à p(t) sin T dx = 2 


then 


It seems reasonable then, that reaction values more nearly equal to those 

resulting from a uniform load will result if Equations (10) and (27) (setting 
= 1) are multiplied by the factor 6/8. Such a procedure was followed in 
this paper. 

Figures 11, 19 and 35 contain plots of the initial elastic reaction expres- 
sion wherein a comparison is made of the curves plotted from Equation (10) for 
the following conditions: 

(1) n 


(2) n= 1 using multiplication factor of “/8 


1 


(3) three terms (n= 1, 3, 5) 

(4) six terms (n = 1, 3, 5, 7, 9, 11) 

This graphical comparison shows that the use of the first mode multiplied by 
the factor e gives results vhich are more conservative approximations to the 
more exact curves. 

Since values of t /T, are determined graphically from the M /Mç versus t/T, 


plot, another source of error in the first mode reaction curves is the use of 








28 
only the first term (n = 1) in Equation (13). However, the elastic center mo- 
ment varles as = for n equal to odd integers, therefore, the approximation for 
center moment a less in error than the corresponding approximation for reaction 
since elastic reaction varies as =. A comparison of the center moments due to 


n 
the two loadings in Figure 36 gives the following result: 


2 
Му" га) = 0.125 steif 


2 
Po) - SE (.500 - .182) = 0.129 p(t) 2° 


then 





M 
(2) e 9 . 082 


By the foregoing reasoning, more accurate values of tm could be gained 
by multiplying Equation (13) by 0.970. This factor is so close to unity, how- 
ever, that little change in the values of t mn would result. In this study, 
Equation (13) with n = 1 was used as derived without applying the correctíon 
factor. 

h. Discussion of Different Phases of Wall Panel Action 

In the introduction to this paper, the statement was made that the charac- 
teristics of the wall panels would cause the shape of the load transmitted to 
the frame to be much different from that which impinges on the panels themselves. 
A study of the Frame Load Charts (end reaction curves) bears out this statement. 
The wall panels are acted upon by an initial-peak, triangular loading pulse 
but send to the frame by means of the end reactions, a loading which has been 


radically altered in shape. 
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The unit load which reaches the supporting members rises gradually from 
zero along the initial elastic reaction curve. Depending on the duration of 
the blast pulse and the beam resistance parameter, the beam (wall panel) be- 
comes plastic at a certain time and the reaction then decreases linearly with 
time at a rate depending only on the pulse duration. Then if the panel has 
sufficient elastic resistance, the action again becomes elastic and remains so 
throughout the remainder of the load pulse. 

For very short pulse durations, the plastic phase continues to the end of 
the pulse and no second elastic phase occurs. Figures 11, 13, and 15 illustrate 
this behavior. As the pulse durations increase, the second elastic phase be- 
gins to appear starting with panels of high elastic resistance (large M/M.). 
Figures 27 and 29 show how (for long pulse durations) the two elastic reaction 
curves oscillate about the incídent load line. It should be noted here that 
damping will tend to decrease the amplitude of the elastic oscillations about 
the load line. If the effects of damping could have been included in the de- 
rivations of the elastic reaction expressions, Figures 27 and 29 would proba- 
bly indicate a gradual decrease with time ín the amplitude of the vibrations of 
elastic reaction about the external load line. Thus, in the second elastic 
phase for low beam resistances, the elastic reaction would become coincident 
with the load line as t approaches EA: 

Theoretically, if a wall panel had infinite elastíc resistance, its end re- 
action would follow the initial elastic curve and the frame loads would rise 
to a maximum. Plastic action takes place however, since the panels have 
finite elastic resistance, and prevents the frame loads from building up to 
these peak values. The formation of the center plastic hinge, thus effectively 
cuts off the top of the elastic reaction curve and consequently permits smaller 


loads to reach the frame. 
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The occurrence of negative values of reaction in the elastic oscillations 
about the load line may be highly significant under certain conditions. If 
the fasteners used to attach the wall panels to the supporting frame are not 
designed to withstand tension loads, it is very possible that such fasteners 
will fail at the time the wall panel end reaction first becomes negative. 
Such failure would alter the frame loads in the case of load pulses of suf- 
ficient duration at, > 0.8 in the frame load charts), limiting the loads 


transmitted by the panels to the supporting members to a single positive pulse. 


B. Tables and Graphs of Results 

Table 1 is presented for convenient reference. For given values of pulse 
duration (Биҝ), plots are also made of the beam resistance parameter (Nol Re ) 
versus both the duration of the initial elastic phase Lol fi and the duration 
of the plastic phase (t/T, - t/T,). (Figures 37 and 39) Similar graphs are 
drawn with M/M. as the fixed parameter. (Figures 38 and 40) 

Figure 37 shows that the duration of the initial elastic phase increases 
almost linearly with increasing values of M/M. for fixed load durations. 
Figure 38 demonstrates that for a given beam resistance, the duration of the 
load pulse has little effect on the length of time which the panel remains 
elastic except in the case of the shorter load pulses. 

As might be expected, there is a very great increase in the duration of the 
Plastic phase (for given pulse duration) with a decrease in wall panel resist- 
ance. (Figure 39) The increase is most pronounced for long load pulse dura- 
tions. 

Figure 40 shows that the effect of load duration is largest for low values 
of panel resistance, and that, in general, the plastic phase duration decreases 


as the duration of the load decreases (with given panel resistance). 
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C. Frame Loads After Termination of Load Pulse 
For time beyond the end of the pulse, the derivation of the elastic reaction 
expression is similar in procedure to that used in Section II in finding the 
limiting reaction for short pulses and constant load impulse. Equation (26) 


would be used to find the displacement and velocity at t = t. in order that 


1 
An and B, could be evaluated; the expression for w would then be found where 


and X. still equals sin nıx/Z. Reaction could be found as before by evaluating 


the shear at x = O. 


VI. SUMMARY AND CONCLUSIONS 
A. Applications of Frame Load Charts 

A method has been found for determining the load-time relation for blast 
load transmitted by the wall panels of a structure to the supporting members. 
The frame load charts may be used for three possible types of wall panel be- 
havior; a purely elastic action, an elasto-plastic action, and a brittle action. 

For a given pulse duration, and purely elastic behavior, the initial elas- 
tic phase reaction curves apply throughout the duration of the pulse. This 
assumes, of course, that the wall panel has sufficient elastic resistance (high 
M /M; ) to behave in this manner. 

A more usual wall panel behavior is the elasto-plastic type. Here the elas- 
tic resistance is not high enough to confine the action to a purely elastic one. 
Given the pulse duration and beam resistance parameter, the frame load charts 
used in conjunction with the dynamic center moment curves indicate the loads 


which are transmitted by the wall panels to the frame. 
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Wall panels constructed of brittle material will behave elastically up to 
the ultimate moment carrying capacity of the panel (My) and then fracture sud- 
denly. The time of fracture could be determineá from the dynamic center moment 
curves by setting Mo = My over a range of values of Mir The frame load charts 
vould indicate a rise in load along the initial elastic curve up to the point 
Of fracture at vhich time the reaction vould drop to zero. 

In order to use (he frame load charts for all three types of vall panel be- 
havior, the following information must be available: 

a. Wall panel span and a curve of side-on overpressure versus time 


from which Pa? t) R., and Ms: can be sica T) 


s? 
b. Dimensions and properties of vall panel in order that T; Mp 

and ку may be calculated or estimated. 
Entering the charts with t,/T, and M/M, or My Ms, the magnitude and time 


distribution of the uniform frame load can be determined, 


B. CONCLUSIONS 

The frame load charts presented in this study represent an attempt to eval- 
uate the loads received by the supporting members from the building covering 
when the walls are subjected to a blast pressure wave. Rough assumptions had 
to be made in some cases in order to obtain fairly simple expressions from which 
the load charts could be constructed. Also, since the elastic reaction expres- 
sions involved the sum of an infinite series of terms, they could not be used 
in their exact form. Therefore, the first mode approximations to the exact 
elastic reaction expressions were modified by a compensating factor of SES and 


then used to construct the frame load charts. It ls felt that this step, 


om M. Newmark, "An Engineering Approach to Blast Resistant Design”, 
Proceedings ASCE, Vol. 79, Separate No. 306, October 1953. 
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although not theoretically exact, provides good engineering approximations to 
the actual blast loads which a frame would experience. Also, keeping in mind 
the additional condition that uncertainties in loading and in the structural 
properties of the wall panels will always exist, the author feels that the 
frame loads indicated by the charts are sufficiently accurate to be used in the 
design of building frames. 

The basic relations for panel end reaction in all phases of elastic and 
plastic action which occur during the load pulse have been derived. For the 
values of pulse duration and beam resistance parameters considered, the load 
charts are camplete through the initial elastic and plastic phases of wall 
panel action, and representative curves for the second elastic phase have been 
drawn. ‘Thus the main purpose of the study -- that of determining the shape of 
the loads transmitted to the frame -- has been accomplished. Additional work 
is required to finish the second elastic phase and to extend the study to time 


beyond the end of the pulse using the methods outlined in this paper. 
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Figure 3— Idealized Elasto-Plastic Moment-Curvature Relation 
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Figure 4—Typical Beam with Loading 
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Figure Y - Free Body Diagram of Plastic Beam 





Figure 8 - Time Origins for Elastic and Plastic Phases 
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Figure 9 - Craphical Representation of Actual and Approximate Variation 
of End Reaction with Center Deflection 
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Figure 10 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versus Time Parameter (t,/T, = 0,1) 
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Figure 11 - Gs of Dynamic End Reaction to Maximum Static End Reaction Versus Time Parameter 
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Figure 12 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versus Time Parameter (t,/T, = 0.2) 
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Figure 13 - ə of Dynamic End Reaction to Maximum Static End Reaction Versus Time Parameter 
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Figure 14 - Ratio of Dynemic Center Moment to Maximum 
Static Center Moment Versus Time Parameter 


(t,/T, = 0.3) 
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Figure 16 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versus Time Parameter (t,/T, = 0.4) 
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Figure 18 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (t,/T, = 0.5) 
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Figure 21 - Ratio of Dynamic End Reaction to Maximum Static End Reaction Versus 
Time Parameter 577 = 0,6 
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Figure 22 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (t,/T, = 0.8) 
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Figure 24 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (6/7, = 1.0) 
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Figure 26 - eye of Dynamic Center Moment to Maximum Static Center Moment Versus Time Parameter 
(t,/T, = 5.0) 
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Figure 28 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (t,/T, = 10.0) 
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Figure 30 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 


Time Parameter (t,/T, = 20.0) 
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Figure 32 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (t,/T, = 50.0) 
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Figure 34 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter (t,/T, = 100.0) 
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Figure 36 - Comparison of Actual and Equivalent Beam Loading 
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Figure 37 - Plot of Duration of Initial Elastic Phase (6/7) Versus Beam Resistance Paremeter (M,/M.) 


for Fixed Values of Pulse Duration (t,/T,) 
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Figure 38 - Plot of Duration of Initial Elastic Phase (t 


T ) Versus Pulse Duration (t,/T,) for 
Fixed Values of Beam Resistance Parameter ( Ms 
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